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ABSTRACT 

A simple and inexpensive way of achieving synchronization in a 
broadband packet network U to ate independent docks at each node 
in the network with the sane mean frequency, and adjust the 
destination clock frequency on the basis of monitoring the destination 
buffer level, which is a linear function of the frequency difference 
between the source and destination clocks. However, due to the 
inherently stochastic nature of packet transport, a periodic packet 
stream of real-time traffic - such as voice, video or a circuit-switched 
signal • tend to become quasi-periodic, causing random fluctuations 
in the destination buffer level. This network-induced jitter makes it 
difficult to extract source clock timing information from the arriving 
packet stream. 

In this paper we illustrate the phenomenon of packet jitter via a 
simple example, and propose methods by which the destination dock 
frequency can be adjusted, so that asymptotically it follows the 
source dock frequency, while minimizing the effects of packet jitter 
on the outgoing packet stream. We provide a comprehensive 
asymptotic analysis of the estimation and control algorithms, and also 
discuss the implementation aspects of these algorithms. 

L, INTRODUCTION 

With the deployment of optical fiber and high speed, 
distributed switching into the public network, broadband packet 
techniques are emerging to be the most viable alternative for 
transporting real- and aonreai-dme signals from diverse traffic 
sources, which range from low speed data to high speed motion 
pictures. In a broadband packet network, information bits from ail 
user services as well is from the existing digital circuit-switched 
traffic are packed into slotted packets, which are transported in a 
truly integrated manner at a network backbone rate, e.g., 150 
Megabits per second (Mbps). Asynchronous Time-Division (ATD) 
Network [1] and Dynamic Time-Division Multiplexed (DTDM) 
Network [2] are two prominent examples of a broadband packet 
network. 

Real-time services generating a periodic packet stream may 
require stringent time transparency, which may call for 
synchronization of the' broadband packet network. Providing 
synchronization through a common muter dock is not a cost- 
effective option, although it may make the transport of a periodic 
packet stream much easier. Alternatively, the arriving periodic 
packet stream at the destination node may be used to provide the 
source dock timing information, and the destination dock may be 
adjusted to follow the source dock frequency by monitoring the 
packet filling level function at the destination buffer, so as to 
reconstruct the original traffic without any loss of information due to 
the frequency difference between the source sad the destination 
docks. Unfortunately, the arriving packet stream at the destination 
buffer is no longer periodic due to random delays introduced by the 
■ network on each packet of the periodic packet stream during transit 
The aperiodidty of the packet arrival at the destination buffer 
appears as a random fluctuation in the destination buffer level, which 
cannot be directly used aa a measure of the frequency difference 
between the source and destination docks. Furthermore, this 
fluctuation may be correlated from packet to packet aa well as 
aonstanonary; its statistics nay vary whenever there would be a 
change in the multiplexing and switching pattern due to addition or 
deletion of calls. 

The central questions in the design of synchronization and dock 
recovery schemes are: (i) how to separate the feeble signal (the small 
frequency difference between the source and the destination clocks) 
from the noise (the packet jitter) with a large magnitude to obtain its 
estimate; and (U) how to control the destination dock frequency 
using the frequency difference estimate. (The terms packet jitter and 
noise are used interchangeably in this paper.) The problem of 
synchronization and dock recovery in broadband packet networks 



has been mentioned in [3] in rather general terms, with no specific 
design details. The reference [4] is the first to hsve attempted to 
address these issues la a quantitative manner. But their assumptions 
are quite restrictive. Moreover, the deadbeat control law used to 
adjust the frequency difference in one time step will accentuate the 
residual estimation error, possibly leading to instability. 

Our approach in this paper is based on the following main 
assumptions: the source dock is generally stable with an essentially 
constant frequency; there is no packet loss during the transit, so that 
the packet stream provides a true, albeit noisy, measure of the source 
dock frequency; packet sequency is unimportant; a bound on the 
packet jitter is known a priori (for the design of the destination 
buffer size); jitter correlation is not considered la the desiga. 

We discuss the problem of dock synchronization in Section 2. 
We first illustrate the phenomenon of packet jitter via a simple 
example of a first-come first-served (FCFS) multiplexer. We develop 
a mathematical model of the destination buffer which describes the 
dynamics of its state, the number of packets in the buffer, as a 
function of the frequency difference between the source and 
destination docks and the random packet jitter. 

Section 3 describes the dock adjustment algorithms, where we 
propose two control algorithms - (1) a first order system, in which 
the current value of the destination dock frequency Is obtained by 
adding its value to the weighted value of the frequency difference 
estimate at the previous instant of time; and (U) a second order 
system, in which the current vtlua of the destination dock frequency 
is derived from a linear combination of the values of the destination 
dock frequency aad the frequency difference estimates at the 
previous two instants of time. 

Section 4 contains two increasingly complex estimation 
algorithms for estimating the frequency difference. la the simplest 
suboptimal algorithm, we use no noise statistics. Therefore, this 
algorithm applies to all situations • packet jitter correlated or not. 
The next level of complexity is added when we explidtly use the 
signal and noise variances in the optimal K aim an filter. Since the 
noise variance is not kaowa a priori, we describe a method to 
process the observation set of the destination buffer state to obtain 
the estimate of the noise variance, too. The latter algorithm is 
optimal with respect to aa uacorrelated packet jitter. 

Ia Section 5, we discuss the implementation aspects of these 
estimation and control algorithms ia a broadband packet network 
environment. Finally, Section 6 contains some conduding remarks. 

2. PROBLEM OP CLOCX SYNCHRONIZATION 

Figure 1 shows a single FCFS packet multiplexer with three 
packet streams with rates 1.0, 2.0 and 4.0 packet per second (pps). 
The packets from these streams compete for the empty slots on the 
high speed line with a rate of 8.0 pps. Ia the sequel we consider the 
4.0 pps rata packet stream to illustrate the packet jitter phenomenon. 

Let f\ be the constant and unknown source clock frequency in 
pps. Let ♦(«). «*0, denote the destination buffer level function. 
The reeding clock frequency at the destination is /j(e) ppe. /,(*) 
changes with « according to the adjustment laws described in the next 
Section. Let A/(«)-/|-/a(«) be the difference between the source 
aad destination dock frequencies. When the packets ia the input line 
arrive with stochastic delays, the input packet stream represents a 
noisy value of the source dock frequency, sad ♦(*) evolves in a 
random manner depending on a/(<) sad the packet jitter 
characteristics. The dynamics of the state of the destination buffer is 
obtained by the conservation law, aad is given by 

*(* ♦ 1) - *(*) + r(«)(A -/i(«)) + «•(*) (2.1 ) 

or 

&«.(*) - <fc(* f 1) - *(*) - r(*)A/ (*) f d («) (2.2) 
where r(*)-l//i(*) is the sampling period and d(t) represents the 
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packet jitter at cth instant. Nonl^^ty in the equation arising from 
tan convenient choice of the sampling period is ignored. (2.1) or 
(2.2) is the state equation of the buffer, based on which all designs in 
this paper are developed. 

Figures 2a, 2b, and 2c, respectively, show the plots of $(*) •» • 
function of * for the 4.0 pps rate packet stream when <i(c)-0 
A/(c)-0, when *(c)#0 and A/(c)-0, and when <*(*)"0 and 
A/(c)*0. The randomness in 6(c) due to contention at the 
multiplexer is obvious from these plots. (Compare Figures 2a and 2b 
at the sampling instants.) 

For a synchronous broadband packet network, synchronized 
through a common master clock. f\ is exactly known and A/(c)-0. 
6(c) in this case fluctuates from e to c + 1 only due to the jitter term. 
Thus, A6(c)-rf(c). (See figure 2b, for example.) 6(c) is a random 
variable. Let max |6(c)|-8, »0 for all sample functions *(*). If 

the first packet from the packet stream is removed from the 
destination buffer after a fixed 8*>8 (allowing for the values of 6(c) 
during the sampling interval) number of packets have accumulated in 
the buffer, then no buffer underflow will occur, and the output 
packet stream from the destination buffer will be uniformly spaced 
with a time period of T-Vf\. 8 is obviously related to the 
maximum aperiodldty introduced into the source packet stream. 

For a plesiochronous broadband packet network, the source and 
destination clock frequencies have the same nominal value but they 
differ by a small amount, which is unknown and whose typical value 
is lew than a few tenth of ppm. Some preliminary studies of 
multiplexer and switch design indicate that the values of </(*) may 
vary from 0 packet to a few packets from one instant to another. 
Thus, the change in 6(«) is mainly due to the random jitter d (t). 
However, since £d(c)-0 (£• is the expectation operator.), we expect 
that the contribution of the A/ (c) term to the buffer level function in 
y(*)T(() time periods may become quite significant, where >(•) 
denotes the length of the interval, which is large and time-varying. 
We propose the following two time-scale model • a faat time 
tfqnrnrrt , A*0, when we collect observation data 6(c), and a slow 
time sequence, m aeO, on which we compute the estimates of the 
frequency difference, denoted by A/(m), from the set 6(*). c*0, 
and adjust /j(m) using this estimated frequency difference. The slow 
time scale interval [m, m + 1) is made time-varying, i.e., Oases/ (m) 
for the above interval. Note that there are a total of j (m)+ 1 samples 
of 6(«) in each large interval, and that the time length of the interval 
[*, a + 11, r(t)-U/j(«). is constant in [m, m + 1) with /»(*)-/»(«•)• 
Obviously, a better estimate of A/(c) would result from a larger 
value of /(•). though this would require a larger destination buffer 
and would cause large delays. We describe later in Section 4 how 
the initial interval, J (0), is chosen and how j(m) is varied with time. 

The most difficult part in the design and analysis of the clock 
recovery schemes is the modeling of the noise process, d{k). Is 
expected value, of course, is zero by the fact that if the frequency 
difference is zero, then the average accumulation of packets in the 
destination buffer must be zero, assuming no packet loss through the 
network. The other statistics, in general, are not known, however. 
d(k) is expected to be a correlated sequence. The correlation is 
dependent on several factors • loading of the network; relative rates 
of the input traffic being multiplexed in comparison to the 
multiplexed line rate; service disciplines at the multiplexers and 
switches, e.g., FCFS or train scheduling; type of connection, e.g., 
virtual circuit or datagram. These factors will Influence the statistics 
of the packet jitter in a complex manner. However, in general the 
correlation of the packet jitter will result from a repeating pattern of 
jitter every few packets in the stream. The period of the pattern is 
large if the packet stream traverses through many stages of 
multiplexing and switching, thereby reducing the jitter correlation. 
The plot of Figure 2b illustrates this repeating pattern of packet jitter 
every 4 packets for the example considered above. 

For simplicity, in this paper we ignore the correlation of the 
process 4(c). We develop an ad hoc design (simple averaging) 
based on no statistics of d(k), and an optimal design (Kalman filter) 
based on d(k) being a white noise sequence. 

3. CLOCK ADJUSTMENT ALGORITHMS 

Since A/(m) is not directly available, we process the 
observation set 6(c), 0*ksj(m) in some manner to filter out the 
effects of the noise and obtain the frequency difference estimate, 
A/(m), for each interval [m, m + 1). We use the estimate A/(m) in 



the control .wXfrithm to compute A (m +1). Let A/(m) denote the 
frequency difference estimation error, i.e., A/(m) - a/(m)- A/(m). 
The estimation algorithms described later are such that the estimate 
A/(m) is unbiased. I.e., £A/(m)-£A/(/ii). Further, if d(k) is a 
white noise sequence, then the resulting A/(m) is also a white noise 
sequence, since it is a linear function of d{k) for all c t [m, m + 1). 
For the purpose of analysis, we assume that. Figure 3 shows the 
signal flow diagram for the estimation and control algorithms, along 
with the two time scales. 

Control Law It The first and simple control taw is to add a 
weighted value of the error-prone A/(m) to /j(m) at m to obtain 
/i(m + l) at m + 1. i.e.. 



/l(* + 1) -/,(/») + aA/(m) 



(3.1) 



for some 0<a<l. With A/(m) and A/(m) defined at above, the 
closed loop system is given by 



/t(m + l) - (l-a)/,(m)+a(/, + A/(m)) 



(3.2) 



To compute the mean of /»(«), we take expectation on both sides of 
the equation (3.2). Since £/i-/i and £A/(m)-0, this results in 



£/i(m + l) - (l-a)£/,(m) + a/, 



(3J) 



The asymptotic value of the mean £/i(») is obtained by taking the 
limit as m -<», which can be shown to yield £/»(») » f\. 

The dynamics of the variance «(m) - E (Ji(m)- Ef^m)) 1 
satisfies the following equation which is derived from (3.2) 

w(m + l) - (l-a)V*) + es 2 <r(i*s) (3.4) 

where o(m)-Bbf*(m) is the variance of the frequency difference 
estimation error. We have assumed that A/(m) is a white noise 
sequence, and thus is uncorrelated with /i(m)-Efi(m) since A/(m) 
is a function of the future values of the packet jitter, which do not 
influence / 2 (m). Since (l-a)<l, the asymptotic variance (as «-•) 

of fi(m) can be shown to be »(») - ■— — »(»). »(•) is taa 

• 2™oi 
asymptotic variance of A/(m), and is related to »(»), the asymptotic 
variance of the packet jitter, d(k), through the estimation algorithms 
described later. The residual variance «(=•) is smaller than the 
asymptotic variance of the frequency difference estimation error for 
all 0<a<l. There is a trade-off between how fast the dynamic 
response should be and how much variance reduction we can obtain. 
A larger a means faster response ((1-a) is closer to zero), but 
higher residual variance, and vice versa. Also, a larger a implies a 
smaller destination buffer. 

Control Law Hi With the control law (3.1), there is not enough 
freedom to control the speed of tracking as well aa the magnitude of 
the residual variance, •»(»), independently. The following control 
algorithm, in which the value of /j(m + 1) is obtained by adding the 
weighted values of fi(m), fi(m - 1), A/(m) and A/(m - 1), provides 
more freedom in controlling the rate of convergence and the 
asymptotic variance. The weightings on A(m) and /j(m - 1) should 
be such that the open loop transfer function should contain (l-x~')* 
as a factor in the denominator (so that there is no steady state error) 
(See [5]). With this, the new control algorithm is given by 

/i(m + l) - (l+8)/ 2 (m)-8/,(»-l) + oA/(m) + eA/(iii-l)(J_5) 

where 8, a and B are free parameters. The dosed loop system can 
be obtained by inserting A/(m)-/ ( -/j(m)+A/(m) into (3.5). 

We will not discuss the complete design procedure here. (See 
Section S.2 of [5] for the design and pp 31 of [6] for choosing the 
optimum parameter values.) It can easily be shown that £/ 2 (oo)-/,. 
The asymptotic variance, «(<■), can be computed from (3.3), and is 
given by 

- M r«» + e»)(l + B + 8) + 2aBfl+8-q) 

v ' l-tt+S-a^+O + BKl+O + Bj + O + SKl+B-cO-te + B) 2 ) 
By making a suitable choice of the two poles of the dosed loop 
system to control the rate of convergence and to minimize •»(•), it 
can be shown that ■»(«) is a monotonic and symmetric function of a 
and 0, and the two poles. Thus, with the poles fixed at the chosen 
values, for »(•) to be minimum, «-fl. The parameter a is then 
determined from the values of the two poles. 

4. FREQUENCY ESTIMATION ALGORITHMS 
The change s in A6(c) from c to c + 1 in each interval [m, m + 1) 
* "i" Is* forwudshift operator, i.e., x»(t)"x(t+l). 




ii mainly due to the noise d (*) than the frequency difference. 

We exploit this separation property in obtaining the estimate of the 
frequency difference from the noisy observation of the rate of change 
of the buffer level function. 

We describe below two increasingly complex algorithms for 
computing A/(m). The first algorithm it the suboptimal simple 
averaging scheme and assumes nothing regarding the jitter process 
except that its ex p ec t ed value is xero. Thus, it is applicable to all 
systems regardless of whether d(Jk) is correlated. The second 
algorithm is the Kalman filter (more accurately, Wiener filter) which 
assumes that rffjt) is a white noise sequence with zero mean and 
finite, but unknown, variance. Since the latter algorithm is optimal 
with respect to a white noise process, it may yield a poor estimate in 
comparison to the simple scheme in case rf(a) is highly correlated. 
While, if <<(*) is white, the Kalman filter will definitely provide a 
better estimate than the simple scheme. Both algorithms are 
developed assuming that control law I is used. These can easily be 
modified for the control law II; we do not discuss it here though. 

Now, if we subtract both sides of (3.1) from /i, then the 
resulting equation gives the dynamic evolution of A/(m) 

A/f>)-(l-«)A/(m-l)-aA/(m-l) «D 

A/(m) responds to its Initial condition A/(0) and the frequency 
difference estimation error 4/(m - 1). The expected value of the 
frequency difference satisfies 

£A/{» - (l-<s)£A/(m) (4.2) 

since £A/f»-0. The initial condition, £1/(0), is known from the 
prior knowledge of the frequency difference. Since A/(m) (and also 
/}(m)) is uncorrelated with the future values of d(k) for 
t • ((it, m + 1) of which A/(m) is a linear function, the variance of 
A/(m) is obtained from (4.1) as 

(l-a) 2 w(m-l) + a 2 o-(m-l) 



«(») 



(4J) 



with the initial variance «(0) known, and o>(m) given. This is the 
same a* (3.4) because w(m)«£(A/(m)-£A/(m)) 2 
£(/j(m)-£/i(m)) 2 , ,inco £/i"/|. Whichever estimation algorithm 
we use in the interval [m, m + 1), the initial mean and variance of 
d/(m) for that interval are given by the above two equations (4.2) 
and (4.3). These quantities are explicitly used in the optimal Kalman 
filtering algorithm described later. 

Computation of Now, we describe how we choose the 

length of the time interval [m, m + 1), i.e., how many dan points we 
use for the estimation of A/(m) in that interval. If j(m) denotes the 
length of the interval [m, m +1), then we propose that 

i(m+1) " 7J^) Km) (4 4) 

with J(0) as the length of the initial period. It can be seen that the 
size of the interval increases at the same rate at which A/(m) decays. 
The rationale behind using this scheme is that as the frequency 
difference decreases with increasing time, the signal to noise ratio 
decreases. Hence, we should have an increasingly larger interval for 
a better estimate of the decreasing A/t>0 with m. 

Both algorithms provide an estimate of A/(m) by processing the 
finite data set A?(«), 0sts/(n) in the interval [m, m + 1). For 
convenience, the buffer state equation (2.2) is rewritten here as 



A*(<) 



(4.3) 



Estimation Algorithm It An ad hoc scheme to estimating A/(m) 
may be just to take the sample mean of all the available observations 
Ad>(-)/i(m). This Is motivated by the idea that the zero mean noise 
samples will cancel out by averaging it. There are a total of y(m) 
samples in the interval [m, m + 1), and thus, 

1 K-o-i 

— 2 *♦(*)/»(»•) <*■«) 



AA») 



•/(») At 



which after simplification gives 



A/(m) - 



Mm) 
Jim) 



(♦"(« + l)-*(m)) 



(4.7) 



where $"(m + 1) is the value of ♦(•) just before the control action is 
taken at the instant m+1. For this algorithm, there is no quantitative 
way to determine how good the estimate is at the end of each 
estimation period, since the noise statistics is not known. However, it 



is possiblef^Fdraw tome conclusions regarding the asymptotic 
behavior of the estimator, if d (a) is assumed to be white. 

„ Mm) >(»>-i 
Consider the error sequence, A/(m) - ■ T d(k), for ail 

m*0. £A/O)-0 since / 2 (m) is independent of </(*) for ail a in 
[m, m + 1) and £d(*)-0. By the same reasoning, the variance of 
estimation error, a(m)-£A/ (m), is given by 

•(«) " ^/"l^ Z '»(*)■ wherc "C*)-" 2 **) is the variance of 
J («•) t-o 

d(k). If the noise variance is bounded, i.e., v(/)sv, for all iaeO and 
for some finite positive constant 7, the estimation error variance is, 
£/f(m)v 

then, bounded by <r(m) * . Now, as j(m)-» with m-», 

J\ m ) 

<r(m)-0, since £/}(«) is a finite quantity. Thus, the frequency 
difference estimation error variance asymptotically goea to zero. 
Asymptotic results are good for analysis. In practice, j(m) is settled 
at some constant value for mxm u where m t is fixed on the basis of 
some criterion discussed in Section 5. »(«) in this case is a finite but 
small quantity. 

Estimation Algorithm tit Further reduction in the variance of 
the estimation error sequence, <r(m), m xO is possible by making use 
of the variances of A/(m) and if (a). The variance, «(m), of A/(m) 
for each m is obtained by the recursive dynamic equation (4.3), since 
ir(0) is assumed to be known and o(m - 1) is known from the 
previous estimation interval [«i -1, m). The noise is assumed to be 
white, with unknown variance. We assume that the noise remains 
stationary in each interval [m, m+1), although, as discussed in the 
Appendix, it is not necessary. The Appendix contains a mechanism 
by which the variance of the noise is estimated in real time using the 
same set of data Ad>(-)- 

Qiuui-itationary Soita Proeou: Let v(m) be the variance of 
d(k) in [m, m + l). Let v(m) be its estimate. The algorithm for 
computing A/(m) is derived in the Appendix, and is given by 

jdsl i 



A/(m) - 



where p(m) 



i+; («•)*(«•) 



(*.-(» + l)-*(m))/,(»i) + 



l+>(»)p(«) 
<t"(m + 1) is as defined above. £A/(m) 



£/j("»)v(i") 

is known from (4.2), which is derived from the known quantity, 
£A/(0). The variance of the error sequence, o>(m), is shown to be 

n(m) 



o-O) - 



(4.8) 



1 + J{m)p{m) 

where «(m) is obtained as explained above via (4.3). v(m) is 
obtained via the methods described in the Appendix (Equation (A.4) 
or (A.5)). This value of o(m) is used in (4.3) to generate «(oi +1). 
The value of EM m ) e>a °* approximated with /j(m) (or with some 
known quantity) (see Section S below) in the above algorithm. 

Approximation with Stationary Soil* Proettt: If the noise 
variance is assumed to be constant, say, v, for all mKO, then the 
algorithm can be simplified as follows. Further, approximate 
£/i(« )"/i(0) (see Section 5 below), since A/(0) is a small quantity. 
Now, p(m)- . In this case, a dynamic equation for p(m) is 

/z(0)7 

obtained from its definition and the equations (4.3) for v(ffi) and 
(4.8) for a(m). 

M^i)-(i-»)'p(")^' lt / ( l;' ( .) «» 

with p (0) - known. Since » is not known, we propose to use 

/!(0)v 

the variance estimation algorithm, (A.4) or (A.5), to estimate the 
constant variance, v, in the interval [0,1), which will be used during 
the other intervals as well. The equation (4.9) is shown to be stable 
for all bounded v(m) and £/j(") in the Appendix. 

It can now be seen that for all bounded p(m) and v(m), 
a(fn)-0 as j(m)-<*> with m-«. Hence, this algorithm, too, has the 
same desired asymptotic behavior as the previous algorithm. It 
should be emphasized that if the sequence d(k) is highly correlated, 
this algorithm may give poor performance. 

S. IMPLEMENTATION OF ALGORITHMS 

The clock frequency is usually expressed as its average value in 
bits per second, along with its tolerance in standard deviation in 
ppm. We will use the packets per second unit for both quantities 
assuming a fixed length packet. Let the source clock frequency be 




/|i/±«„ where / U the nomw^pfequency and t, ii it* standard 
deviation. Likewise,- let the initial destination clock frequency be 
/i(0) -f±U- ThU suggests that £/i-E/ 2 (0)-/. and the variance of 
/, u «? and that of / 2 (0) U ti. Therefore, £A/(0)-0 and 
«(0)-tJ+«J- 

Now, the nest question is how to choose >(0)? J(0) i» 
determined from the prior knowledge of A/(0) (or «(0)), the 
maximum value of the noise magnitude and the allowable buffer size, 
so that even if the system is running for J(0) units of time in the 
presence of noise, there is no buffer under- or overflow, and we get 
a reasonable estimate of A/(0) at the end of the first interval. In 
practice, after a few iterations A/(m) will become very small. At 
this time we can stop changing j (m) and use a constant value. We 
propose that for some m -mi, whenever either |A/(mi)|<ii, or 

«(m.)<«l. or 'j 1 "'? <«i, then set j(m)-j(m + 1) for all m*m,, 

where t lt < 2 and tj are design parameters. 

Rtutttng of Algorithms: Resetting may be necessitated whenever 
docks at a node are switched or replaced. We propose the following 
mechanism. Any time the value of the sample function |<t(*) |*p, we 
reset the algorithm with the same set of initial conditions. The 
assumption is that the clocks have the same kind of frequency 
tolerance, sp is the upper and lower limits on the buffer size, and 
can be determined on the basis of the maximum expected deviations 
in ♦(*) due to the packet jitter, the expected values of 
A/(mi), /aO»t), and the changed frequency difference. As long as 
we can choose all the parameters conservatively to avoid buffer over- 
or underflow, a false resetting of the algorithms makes little 
difference in the steady state. 

In the above algorithms, the time intervals [m, m +• 1), m »0 
are determined a priori, and the destination buffer levels are read at 
time instants, m aO, the values of which may contain a fraction of 
packets. Alternatively, it is also possible to fix the destination buffer 
level function at discrete packet levels, e.g., half, three quarter etc, 
and to measure the time intervals for ♦(«) to attain these 
predetermined labeled values. There are two limitations of this 
technique: (i) 6(e) may occasionally cross the labeled points due to a 
large but unexpected jitter. If for some reason these crossings occur 
too frequently, the estimation of the frequency difference may be 
quite poor due to small number of observations, which will result in 
a larger residual jitter. (11) By observing 6(e) at a coarse level, we 
loose information about the statistics of the high frequency jitter 
component buried in it. Thus, the variance of the jitter cannot be 
estimated on-line from this observation set, and the Kalman 
algorithm cannot be used. 

If i(t) u highly correlated, then the control law a with the 
estimation algorithm I would be a preferred combination. This is 
because the residual effects on / 2 (m) due to the poor estimates 
A/0") can be further filtered out better by the second control law 
than the first one. 

6. CONCLUDING REMARKS 

We have developed clock recovery schemes for the real-time 
traffic in a broadband packet network. The algorithms for the 
adjustment of the clock frequency are based on a mathematical model 
of the destination buffer, and the assumptions made regarding the 
nature of the packet jitter introduced by the network on a periodic 
packet stream generated by the source. We process the tame 
observation set, the destination buffer level function, to obtain the 
estimates of both quantities - the frequency difference between the 
source and destination clocks, and the jitter statistics; we perform 
some sort of averaging to obtain the former, while a sort of 
differencing yields the estimate of the jitter variance. It has been 
analytically shown that the destination clock frequency adjusted 
according to the proposed algorithms follows the source clock 
frequency asymptotically, with a minimum of fluctuation. 

We have provided a systematic design procedure for choosing 
the relevant parameters for implementation, and have alto discussed 
the sensitivity of the performance with respect to these parameters. 
We plan to test the validity and effectiveness of the algorithms via 
simulation by Implementing them on a computer model of a single 
stage FCFS multiplexer for different source traffic rates, which we 
will report in a future paper. 
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APPENDIX 

Kalman Filter. The Kalman filter is the best example of an 
optimal filtering algorithm based on the mean square criterion and a 
white noise sequence. The reference [7] provides a good tutorial 
introduction to the Kalman filter. 

A/0*) and /j(m) are constant in the interval [m, m + 1). Thus, 
for the development of the Kalman filter algorithm thai state and 
observation equations in that interval are given by 

A/t« + l)-AA«) 

Aq(«)/ 2 (») - A/(*)+/ 2 (»)rf(*) 

The Kalman filter equation with the estimation error, A/(*), 
(A/(».)-A/(*)-A/(«))ia 

A/(« + l) - A/W + *(* + l)(A*(« + l>/ 2 (m)-A/(*)) (A.1) 

with 

A/(0) - £A/(m) + i (0)(A*(0)/ 2 (m)-EA/(»)) (AJi. 

as the initial condition, where £A/(m) - (l-a) a £A/(0) is the 
expected value of the initial frequency difference, which is known. 
The gain, «(*). and the estimation error variance, e>(*)-fc"A/ (*), 
are given by 

with <r(0) - £A/*(0) - £(A/(m)-£A/(m)) 2 - «(m) is the initial 
condition, which is also available from the recursive equation (4.3). 
£/f(m)v(a) is the variance of /j(m)<f (i), since / 2 (m) is independent 
of d(k) for ail k * [m, m + l), and £<f (*)-0. 

If the noise is nonstationary, then (A.l), (A.2) and (A.3) 
describe the recursive Kalman filter algorithm with »(*) replaced 
with its estimate, f (*), obtained through the methods described later. 
For simplicity, we assume that the noise is piecewise nonstationary, 
i.e., its variance is constant in each interval [n, m +1). In this case, 
the filter equation can be solved analytically to obtain a simple 
nonrecursive estimator to compute A/(m). 

The analytic solution of the above difference equation for »(«) 



»(*) - 



«t(0) 



1 + 



_2i°L 



£7i(m)S(m) 



with o(0)»w(m). v(m) denotes the estimate of the constant noise 
variance in the interval [m, » + l). Now, with this o(t) and i («) 
from above, the filter equation (A.l) can be solved to yield 

P(m) ,_«,/-» ^ l 



A/(«) 



1+./ (»•)*(«) 



(*-(» + l)-*(m))/ 2 (m) + 



l+>(«)p(m) 



£A/(m 




with p (m)- 2 ^ ■ ♦"(■^^J >» the value of the buffer level 

?/!("»)»('") 

function at the insunt ju»t before m + 1, and w(in), 5(m), £A/(m), 
;(m) and ire all known quantities. For the purpoie of 

computation, we will take £/i(« WiO")- 

On-lUu Computation ofScUt V or lata Estimatt. v(t): Now, we 
develop * method for estimating v(c) from the available information, 
a*(t + l). A*(«) A*.(0). 

Approximate Anatytls: Since A/(«) U small, we can ignore its 
effect on A*(«) on the fast time scale. Therefore, a*(*)-d(«) for 
all t*0. Hence, v(«) - £d 2 («) - EA* 2 («). Since the probability 
density function of A6(c) Is unknown, the variance, »(*), can be 
estimated as follows (»(*) is the estimate of v(t) at «.) 

for all «>*aeO, with 5(0)-A* 2 (0), or some initial guess of 5(0), if 
it is available. 

Exact AaatytU: Since Mm) and, therefore, a/(m) are both 
constant in the interval [m, « + l), successive differencing of the 
observation A6(*) yields 

AA*(*) - A*(« + l)-A(««)-d(* + l)-<i(«) 

The expected value of AA*(<) is zero, and its variance is given by 
£AA* , (*)-eT(*'(* + l)-d (*))*. Since the noUe is uncorrelated. this 
yields 

»(«••■ l)+v(«) - £d 2 (« + l)+£d 2 (*) - £AA* 2 («) 

v (i) in general is time-varying, but if we assume that its value does 
not change significantly from * to * + 1, then v (* + !)«¥(*)• Thus, 



v(*)-£d 2 (*)-y£AA* 2 (*) 



As before, 5(4) can be computed from the following recursive 
equation 

k J_J. flA(> 2 (t) 



« + l 2 



(A J) 



- AA^ 2 (0) u ^ iaitjai condition. The 



with 5(0) 

computation starts at m -0, and must be frozen at time instants /(«•) 
and /(m)-l because in the definition of aa*() two future values of 
♦(•) are needed and the value of A/Cm) changes at the dme instant m 
due to the control action. Here we assume that the effect of this 
control action on ♦(•) is observed only at the next time instant, i.e., 
one fast time-scale dme period later. 

If the noise statistics is such that the variance v(-) is a constant, 
then the recursive variance estimator will come to a steady state 
value quite soon, perhaps in the first slow dme scale interval. At this 
dme, the estimation of the noise variance can stop. Which variance 
estimation algorithm should be used will depend on the relative 
magnitude of a/(0) and the expected packet jitter. 

ProbUm of Stability: There is a potential problem of stability 
due to the time-varying nature of the variance estimator, and <x(m) 
being generated by the nonlinear equation (A.3) and used in the 
computation of «(m + l) via (4.3). 5(m) computed via (A.4) or 
(A.S) is a bounded quantity, i.e., 0<5(»)ssv, for some constant v. 
£/}(«■) is also a bounded quantity, i.e., |£/t(m)|*7i. tot tome 
finite constant/}. Then it is easy to see that, with p (m)- , the 

equation for p (si) is obtained from (4.3) as 

f( " M11(tl - )1+ ^w>(.) );W 

It is not difficult to show that if 0<a<l, then the coefficient of p («) 
in the above equation is less than one for all p(m), bounded or 
unbounded. Hence, p (m)-0 as ;(m)-» with m-«. 
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Rpm 2 : Buffer FUIiaf Laval at a Function of Timt 
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